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We investigate time- dependent properties of a single particle model in which a 
random walker moves on a triangle and is subjected to non-local boundary conditions. 
This model exhibits spontaneous breaking of a symmetry. The reduced size of 
the configuration space (compared to related many-particle models that also show 
spontaneous symmetry breaking) allows us to study the spectrum of the time-evolution 
operator. We break the symmetry explicitly and find a stable phase, and a meta-stable 
phase which vanishes at a spinodal point. At this point, the spectrum of the time 
evolution operator has a gapless and universal band of excitations with a dynamical 
critical exponent z — 1. Surprisingly, the imaginary parts of the eigenvalues Ej(L) 
are equally spaced, following the rule ImEj(L) oc j/L. Away from the spinodal point, 
we find two time scales in the spectrum. These results are related to scaling functions 
for the mean path of the random walker and to first passage times. For the spinodal 
point, we find universal scaling behavior. A simplified version of the model which can 
be handled analytically is also presented. 
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1. Introduction 

Spontaneous symmetry breaking in non-equilibrium statistical mechanics was recently 
observed in several one-dimensional many-particle models [[I], 0, |3], |], [5|]. In this paper we 
consider a single particle random walker model || which exhibits symmetry breaking 
of a Z 2 symmetry. We investigate the spectrum of the time evolution operator, and 
stationary and time-dependent properties. The random walker moves in the two- 
dimensional geometry of a discretized right-angled triangle. In the interior of the triangle 
it may hop locally to neighbouring sites, whereas on the two short sides it may jump 
non-locally to a corner. Unusually for a one particle model the simple model, studied 
here shows spontaneous symmetry breaking. Considering the infinite volume limit, the 
particle stays in one corner of the triangle and no longer reaches the opposite corner 
which is reflected under the symmetry. This yields two equally stable phases. 

Breaking the symmetry explicitly in a soft way, one of the stable phases first becomes 
meta-stable, before one reaches a spinodal point where it becomes unstable. Hence one 
can investigate metastability and spinodal points in non-equilibrium processes with the 
simple random walker model. 

The random walker on a triangle (RWT) model is defined as follows. The walker 
moves on the triangle 

Tjj = {(j, k) | j, k integers > with j + k < L} . 

We define the following processes for the random walker at site (J, k) (see also Ref. 0), 
where p is the probability of a transition in an infinitesimal time step dt: 

U, k) -> (j - 1, k) if j > 1 with p = a(l + h)dt 

(j, k) — > (j, k - 1) if k > 1 with p = a(l - h)dt 

(j, k) —> (j — l,k+ 1) if j > 1 with p = 6(1 + h)dt 

(j, k) -> (j + 1, k - 1) if k > 1 with p = b(l- h)dt 

(l,k) -> (1,L- 1) with p = c(l + h)dt 

(j,l) -> (L- 1,1) with p = c(l - h)dt . (1.1) 

The first four processes are local hops to neighbouring sites on the triangle, the last 
two are non-local jumps to the corners (1, L — 1) and (L — 1, 1) respectively. For h = 
the definitions are invariant under reflection on the line j = k. We fix the unit of time 
choosing c = 1/2. This renders the probability rates dimensionless. 

In the following we also impose the condition a + b = 1/2, which makes the 
RWT model the "zero-temperature" limit of a three-state model |], |J which shows 
spontaneous CP symmetry breaking. (Therefore the reflection symmetry of the RWT 
model will also be denoted as CP symmetry in the following.) It would be interesting 
to study the general case a + 6^1/2as well. In the three-state model, positive and 
negative particles hop among vacancies on a one-dimensional chain. The C operation 
interchanges positive and negative particles, and the P operation interchanges left and 
right on the chain. The investigation of the free energy functional, and of the flip times 
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Figure 1. Simplified phase diagram of the three-state model. In the "low- 
temperature" regime A the symmetry is spontaneously broken for h = where h 
is the strength of a symmetry breaking field. It is separated by a line of spinodal 
points from the "high-temperature" regime B. 



between stable and meta-stable phases, leads to the phase structure given in figure [I] 
||. In the CP-symmetric case, one finds two broken phases J/]] for "low temperatures" 
(in this model the output rates of the particles play the role of a temperature) and a 
disordered phase for "temperatures" above T c . At the critical point (T = T c , h — 0), one 
finds a dynamical critical exponent z — 2, which is connected with the appearance of 
shocks || [10], |H| . Breaking the symmetry explicitly, one finds a "low-temperature" 
regime A with a stable and a metastable phase. This regime is separated by a spinodal 
line from the regime B (Fig.|]). At the spinodal line, a previous investigation suggested 
z = 1 for h 7^ |§. This motivated us to study time-dependent properties, metastability 
and spinodal points in more detail. All these are manifestly associated with non- 



equilibrium models |L2 



Another model which shows spontaneous CP symmetry breaking is the two-state 
(particles and vacancies) asymmetric exclusion model 0, [3], |J. There, the C operation 
interchanges particles and vacancies. In the spontaneously broken stationary state, 
shock profiles (connected with a dynamical critical exponent z = 2) appear. But 
breaking the symmetry in a soft way one does not find a region with a meta-stable 
phase, and hence there are no spinodal points present in the model. 



The time evolution of the RWT model is given by the master equation [13 



±\P>=-H\P> , (1.2) 

where \P> is the vector of the probabilities P(j, k) of finding the random walker at site 
(j, k) . The "hamiltonian" H is the time evolution operator and is determined by the 



processes (|1 . 1| ) . The stationary properties of the system are given by the ground state 
\Po> with H \ Pq> = 0, and the dynamical properties are given by the excitations of H . 

The configuration space of many-particle models grows exponentially with the 
system size in contrast to the configuration space of the RWT model which has dimension 
L(L — l)/2. This polynomial dependence on L enabled us to investigate the spectrum 
of the model. By numerical diagonalization of the hamiltonian, we found an interesting 
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two band structure of the spectrum, which leads to two gaps or time-scales. At the 
spinodal point h c = a/(l — a), one scale vanishes with critical exponent x = 2 and the 
lower part of the spectrum has the form Ej = akj/L. The constants kj are universal, 
and do not depend on the remaining parameter a. This reminds one of the critical 
spectrum of a quantum spin chain related to a conformal field theory. As a surprise, the 
imaginary parts of the constants kj have the form Imkj = jlmki. For the real parts, we 
do not find any regularity. 

The value of h c has already been conjectured in ||. It is also predicted by taking the 
"zero-temperature" limit of the mean-field results for the three-state model discussed 
in ||. There, mean-field theory was found to be compatible with simulations for small 
"temperatures" (output rates). 

Besides the study of the time evolution operator one can, as with other models, 
investigate the free energy functional (FEF) of the RWT model. This functional / is 



where d is an order parameter of the model and Po(d,L) denotes the corresponding 
stationary probability distribution ||. For the RWT model, we take 



as the order parameter. 

Analogously to the "low-temperature" regime of the three-state model, one finds 
that the FEF of the RWT model has two global minima (d = —1,-1-1) in the CP- 
symmetric case (h = 0). In the limit L — > oo, the symmetry is spontaneously broken 
and there are two equally stable phases with d — ±1. For < h < h c , the dynamics 
( |1.1|) favors the movement to lower values of the order parameter and the FEF has a 
global minimum (d = —1) and a local minimum (d = +1). This leads to the existence 
of a stable and a meta-stable phase in the infinite volume limit (L — ► oo), the random 
walker staying near the favored corner F=(l, L—l) and the unfavored corner U=(L— 1, 1) 
respectively. Above h c there exists only one phase, the random walker always staying 
around F. This is also reflected in the spectrum of the hamiltonian ( |1.2| ), since one has 
two ground states for L — > oo and h < h c and only one for h > h c . 

It is also possible to measure the flip time from one phase to the other using Monte 
Carlo simulations. At the spinodal point, one observes a discontinuity in the passage 
time from the corner U to the corner F. At and above h c , the average time taken 
to reach the corner F grows linearly with L. However, there is a discontinuity in the 
proportionality constant as one approaches h c from above. The linearity above h c can be 
explained by a mean- field argument, the linearity at h c cannot. The time taken to reach 
the corner U and the time taken to reach the corner F for h < h c grow exponentially 
with L. 

Using Monte Carlo simulations, one can also investigate the movement of the random 
walker starting from specific initial conditions. Analysing the mean path, one again 



defined as P, |H|: 



h(d) = --logP (d,L) , f(d) = Km f L (d) , 



(1.3) 



d = (j - k)/L 



(1.4) 
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identifies two time scales away from the spinodal point, and at the spinodal point, one 
again finds universal behavior. 

All the results found here for the RWT model are in accordance with the results 
found for the three-state and the asymmetric diffusion models. But the detailed analysis 
of the spectrum of the time evolution operator now suggests a link between the time- 
scales found in the spectrum, the time-scales found in the flip times, and the scales one 
observes in the movement of the random walker. This was not observed before. 

The random walker on a triangle can be simplified to a random walker on a right 
angle (RWRA) model which can be handled analytically while still possessing a meta- 
stable phase. This is done taking b = and keeping only the points (j, k) G 72 with 
j — 1 or k — 1. That is, we consider only the movement along the j- and /c-axes, where 
macroscopic jumps occur. In an infinitesimal time interval dt the following processes 
can occur with a probability p: 



The RWRA model is essentially a one-dimensional (albeit non-local) random walker 
model. Again we limit our investigation to the case c = 1/2 (this fixes the unit of time) 
and a = a + b = 1/2. We checked that taking a different from 1/2 does not change the 
physics of the RWRA model. One is able to calculate exactly a characteristic equation 
for the eigenvalues of the time evolution operator, the free energy functional and the 
flip times. The spinodal point is at h c = 1. Since all rates in fll.5|) have to be positive 
one can study with the RWRA model only the regime below the spinodal point. The 
results are given in appendix C. 

The paper is organized as follows. In the next section, we present our results on the 
spectrum of the hamiltonian of the RWT model. We identify two time scales, one of 
which vanishes at the spinodal point. At this point, the low-lying excitations acquire a 
universal L~ l dependence yielding z = 1. We also show that finite-size scaling relations 
hold. 

In section 3, we compute the FEF of the RWT model and present our Monte Carlo 
results on the passage times from the favored corner F to the unfavored corner U and 
vice versa. Two scales can be measured and related to the depths of the minima of the 
FEF and to gaps in the spectrum of H found in section 2. 

In section 4, we give results on the movement of the random walker that mirror the 
scales found for the spectrum (section 2). At the spinodal point, one finds a non-trivial 
universal behavior that also reflects the dynamical exponent z = 1. Section 5 contains 
our conclusions. 

In appendix A, we give some simple results for the movement of the random walker 
for h above the spinodal point. In appendix B we explain the connection between the 
flip times and the spectra of models with absorbing configurations. In appendix C we 
investigate the random walker on a right angle model. 



U, 1) 

(1,*) 
(l,fc) 



(7-1,1) 
(£-M) 

(M-i) 

(1,L-1) 



if j > 1 with p = a(l + h)dt 
with p = c(l — h)dt 
if k > 1 with p = a(l — h)dt 
with p = c(l + h)dt 



(1.5) 
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Figure 2. Dependence of the spectrum on system size. The data is from numerical 
diagonalization for a = j, h — 0.3. Shown are the real parts of the eigenvalues of H. 
The gaps m\ and extrapolate to finite values as seen in figure |3[ 



2. Spectrum of the hamiltonian 

Using a modified version of the Arnoldi algorithm [13], we have calculated the low- 



lying eigenvalues (up to the 30th level) of the time evolution operator H for L < 220, 
various values of h, and a = 0.001, 0.1, 0.25, 0.4. Apart from the ground state "energy" 
E (L) = 0, the further eigenvalues of the non-hermitian matrix H are positive, or 
complex with positive real parts. Since the hamiltonian matrix has real entries only, 
complex eigenvalues come in conjugated pairs. Typical results for the real part of the 
spectrum (for a = 0.25, h = 0.3) are shown in figure One identifies two bands of 
excitations leading to the gaps (time scales) mi and m2 in the large L limit. Note that 
the excitations correspond to energy gaps, since the ground state energy of H is zero 
for any system size. 

Values for the large L limit were extrapolated from data for finite L using the 
Bulirsch-Stoer algorithm [ IB ]. 



The most important observation is that there exists a critical spinodal point at 
h c = (see figure 3). At this point, the system loses one of its time scales (m-i = 0) 
and acquires a band of excitations with a universal L~ l dependence. 

The results for the spectrum of H can be summarized as follows: 

• first excitation: 

For h < h c the first excited energy level approaches zero as 

Ei oc exp(-tuL) (2.1) 
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Figure 3. Length scale u, and upper and lower gaps mi and m,2 for a = The values 
given are extrapolated from data of numerical diagonalization of H up to L = 220. 



This excitation is real, the higher excitations are complex conjugated pairs. Thus, 
in the limit L — > oo, there are two ground states. Linear combinations of the 
eigenvectors give two phases J7J. For h — they are equally stable, while for h > 
one of them is stable, and the other is meta-stable. The two phases correspond to 
minima of the free energy functional at d = 1 and d = —1 (see section 3). The 
meta-stable phase vanishes at the spinodal point h c . Above h c , there is only one 
ground state for L —>■ oo, i.e. none of the excited levels approaches zero. 

• higher excitations: 

Above Eq and E±, there are two bands of complex eigenvalues of the form 

Ef\L) = mi + kf ) /L + --- (2.2) 

Ef\L)=m 2 + kf ) /L + --- . (2.3) 

The gaps mi and m 2 are real, whereas the constants have imaginary parts, 
such that one has complex conjugated pairs of eigenvalues. One finds that mi is 
proportional to u. For a = 0.25. we find mi = 0.57(4) u. 

For h = one has mi = m 2 and the constants = k^ = k are independent 
of j, i.e. the 1/L term is the same for all excitations in contrast to Eqs.( [2.2| ) and 
(|2.3j ) for h > 0. There are two ground states (compare Eq. fl2.1|) ), corresponding to 
equally stable phases. 

For h > h c the excitations again come in the form of two bands as in Eqs. fl2.2p and 
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k ■ 

Table 1. Excited levels of H at h c : Ej = a-jj. Complex excitations come in conjugated 
pairs. 



level 


Re(^ (1) ) 




Im(A:j 1) )/12.81 


1 


1.7080 








2 


8.4980 


12.7663 


0.997 


4 


10.368 


25.771 


2.012 


6 


11.504 


38.548 


3.009 


8 


12.324 


51.241 


4.000 


10 


12.97 


63.89 


4.988 


12 


13.5 


76.52 


5.994 



lower gap: 

For h — > h c , the lower gap ra\ vanishes, as does u, which is proportional to m\. 
For h < h c , the gap is given by ra\ = lini^oo E 2 (L). Above h c , the first excitation 
does not give a second ground state for L — > oo and one finds a non-zero gap 
mi = limi^oo Ei(L). The gap mi given by 

mi = 0.42(2) (h e - h) x for h < h c 

m x = 0.46(2) (h - h c ) x for h > h c (2.4) 

where the numbers are given for a = 0.25. Note that these equations are valid for 
any h and not only in the vicinity of h c . One finds an exponent x = 2.00(1) for 
h < h c and h > h c . Extrapolated values for the m ; are shown in figure |3|. 
The position of the minimum of the parabola Eq. (|2.4| ) determines the spinodal 
point at h c = with a precision of two digits. 

upper gap: 

The scale 777,2 increases monotonically with h over the entire range < h < 1. One 
finds 

m 2 = 0.48(2)/7 L22(2) + 0.043(2) (2.5) 

for a = 0.25. 

lower band of excitations at h = h c : 

At the spinodal point, the lower band is of the following form: 
E = 

^ fci const. 

1 = a L + + " ' ( ^ 

_ kj const. . _ . , 

hj = a— 1 ^ — h • • • 7 > 2, complex conjugated pairs, 

with the extrapolated values kj given in table |l|. One has universality: apart from 
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the normalization constant a ("speed of light"), the leading terms are the same 
for any critical point h c (a) = Extrapolation of the k^ for a =0.25, 0.4, 0.1 
gives the same result with 5 digits precision. At the spinodal phase transition, the 
dynamical critical exponent |T7[ is z — 1, as was found for the three-state diffusion 
model where, however, this was only seen for the first excitation ||. 
The imaginary parts of the numbers kj follow (see table P 

Im(kj^) = 12.81(1) j . (2.7) 

This equidistant spacing comes as a surprise and we do not know its reason. For 
quantum spin chains related to a conformal field theory one expects regularities 
(conformal towers) for the real parts of the kj. However, we could not detect any 
such structure looking at the real part of the spectrum. This should be a result 
of the reduced state space of the RWT model. Because of the large configuration 
space of the full three-state model, a similar analysis was not possible there. 
For the higher band of excitations at the spinodal point, one again has the form of 
Eq.(P). 

finite-size scaling analysis: 

Finite-size scaling relations, which describe the behavior of finite quantum spin 



chains in the vicinity of an (equilibrium) critical point of the infinite system [18 
apply also to the (non-equilibrium) RWT model at h c . Consider the curves 

F L : LE 1 (h,L) . (2.8) 

The critical point h c can be determined from extrapolation of the crossing points 
Across {L) of the curves Fl and Fl-\. Indeed one finds 

lim Across 

(L) = 0.3333(2) = h c (2.9) 

L^oo 

for a = 0.25. This is the most precise numerical determination of h c . 

Following the finite-size scaling hypothesis, one can also determine the critical 



exponent x in a way different from that of Eq. (|2.4|) ||18|| . From the available data 
up to L = 220, one again finds x = 2.0(2). 

In the following sections, we discuss further properties of the RWT model that can be 
related to the scales of the spectrum and to the universality at the critical point. The 
scale uj can be linked with flip times and the stationary FEF (see next section). In 
section 4, we present dynamic quantities related to the properties of H, specifically to 
the gaps ffij. 



3. Free energy functional and flip times between phases 

We have investigated the behavior of the stationary probabilities Po(j, k) for different h. 
These probabilities were determined from the stationary state as computed by numerical 
diagonalization of the hamiltonian (see section 2). First, we remark that for h = all 
eigenstates are also eigenstates of the CP operation. We find that the eigenvectors 
corresponding to the complex conjugated pairs of eigenvalues Ej of the lower band of 
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d 



Figure 4. The free energy functional /l(c?) for a = j, L = 100 and h = 
0.0,0.1,0.2,0.3333,0.4. The values are shifted by the value The data is 

generated by numerical diagonalization of the hamiltonian. Since the precision is 
limited to 10~ 12 the curves for larger h cannot be fully determined. 



1e+08 



20 40 60 80 100 120 140 160 180 200 



Figure 5. The flip times T s h or t(-^) and Ti ong (L) for a = j and h = 0.1. The fitted 
lines are given by 50exp(0.2L) and 500 exp(0.052 L). 
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excitations ( |2.2| ) are CP + for j even and CP~ for j odd. 

One also finds that, for any h and finite L, the eigenvector \P > with eigenvalue 
= is the only vector that directly corresponds to a probability distribution, i.e. all 
its entries are positive. In the case h = 0, \Pq> has CP + symmetry. 

Hence, for h = 0, the free energy functional f(d), defined in Eq. fll.3p , is symmetric 
under the inflection d — > —d. (From the data, one finds that the free energy functional 
is a well defined limit as was already known for the asymmetric exclusion model and 
the three-state diffusion model [§].) There are minima at d — — 1 and d — 1 and a 
maximum at d — 0. The minima correspond to the two equally stable phases (see 
figure ^|). This is basically the same situation as for the CP-symmetric three-state 
model in the "low-temperature" phase. 

For < h < h c , one has a global minimum at d — — 1 and a further, shallower 
minimum at d = 1 (see figure |J). This reflects the fact that one of the phases becomes 
meta-stable. One can define two scales: 

^large /max jf( 1) 

^small = /max ~ /(I) , (3.1) 

where / max is the maximum value of f(d). We shall relate these scales to the time the 
system needs to flip from one minimum to the other. At h c , the second minimum of / 
(and hence the meta-stable phase) vanishes. This cannot be seen in figure |] but can be 
inferred from the results for smaller lattice sizes. For h > h c , the free energy has only 
one minimum at d — —1. 

Analysis of the eigenvectors of H that become the steady states in the large L limit, 
and the relation of the RWT to the three-state diffusion model, have given the result 
that below h c there are two distinct phases with order parameter values d — 1 and 
d = — 1. Above h c , however, there is only one phase: the walker stays near the favored 
corner F of 71 . This can be verified with Monte Carlo simulations. In a way similar 
to that for three-state diffusion model [p], ||, one can measure the time for the RWT 
model to flip from one phase to the other. One finds that, for large L and h < h c , in 
the neighbourhood of a corner F=(l, L — 1) or U=(L — 1,1), the random walker moves 
through the configurations on a shorter time scale than is needed to get from one corner 
to the other. This allows one to measure the flip times Ti ong from favored to unfavored 
and T s hort from unfavored to favored phase as average first passage times from F to U 
and U to F respectively. 

With c = \ and a + b = |, the unit of time is defined such that, on average, the 
random walker moves one step (|1.1 ) per unit. From the simulations, one finds that, in 



leading order, the flip times are given by 

Ti ong oc exp(/ii ong L) for any h (3.2) 

J exp(/i short L) for h < h c , . 

lshoit ^{L ioih>h c . {6 - 6) 

For h < h c , the time intervals the random walker spends alternately near either corner 
increase exponentially with L. The fraction of time in the meta-stable region decreases 
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exponentially. In this sense, for h < h c , there are a meta-stable and a stable phase. 
Results of the measurement for a = 0.25 and h = 0.1 are presented in figure [5[ 

For h > h c , the proportionality constant in Eq. (pT3|) can be calculated by a simple 
argument, see appendix A: 

T h>h c _ L _ 1 + he , v 

J sho r t - h + a{1 _ h) ~ h + hc L - ^ 

In contrast, the appearance of a linearly growing flip time at h — h c is non-trivial. 
One finds the form 

T short = - - — — + const. , (3.5) 
a a 

where 7 = 0.57(6) is a positive constant and independent of a to a precision of 
about 10%. Up to the normalization a, the short flip time is, in leading order, a 



universal function of L as are the leading terms of the lower spectrum Eq. (|2.6| ). The 
appearance of a universal linear flip time is linked to the universal 1/L spectrum of 
the hamiltonian, compare Eq. (|2.6| ). One can ask how general this universal behavior is. 
Similar investigations for other models remain to be done. 

Note that from Q3.4|) and (|3.5| ) one finds that in leading order 

Jim T*>£ = 2T s t rt • (3.6) 



h—>h c 



Hence there is a discontinuity in the flip time T^ OTt . At the spinodal point, T s ^ c ort has a 
special L dependence, essentially different from its behavior above or below this point. 

The values of z^maii and vi arge in ( |3.1| ) , determined from numerical diagonalization of 
H for L < 140, and measurement of the flip times, up to L = 200 for T s h or t and L = 40 
for Ti on g, yield the following result: 

tshort = ^small for h < tl c 

AHong = ^large for any h . (3.7) 

As an example for h = 0.1 and a = 0.25, we find the values /i s hort = 0.052(3) 
and /iiong = 0.20(2) (see Fig.|) which is to be compared with z/ small = 0.053(2) and 
^large = 0.17(2) (see Fig.f|). For h = and a = 0.25 the equations ( |3.7| ) can be checked 
to three digits; the precision decreases with increasing h to about one digit at h = 0.3. 
(Finite-size corrections to ( |3.1| ), ( |3.2| ) and ( |3.3| ) become more important for higher h due 
to crossover effects to the linear regime for h > h c .) 

Equation ( |3.7| ) is interesting because it relates (steady state) properties of the free 
energy functional f(d) with (dynamical) flips between phases. This is unexpected for 
a non-local model without detailed balance. For the RWRA model (appendix C) an 
analogous link can be made analytically. Also, for the three-state diffusion model, 
similar observations, albeit with less precision, were made in ||. 

One also finds that /i s hort = ^ (recall that E x oc exp(— ojL) is the first excitation of 
the hamiltonian, see Eq. ([2.1|) ). For h = 0.1 and a = 0.25 the value of u from numerical 
diagonalization is 0.048(3). The reason for this connection, together with an analogous 
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s(t)L 




t 



Figure 6. Average distance s(t) from starting point U=(L — 1, 1). 

Data from simulation of the RWT model for a = j, h = 0.3 and system sizes 

L = 100, 200, 400, 600. The scaling variable is t. 

expression for /ii ong , is given in appendix B, where we consider modified models with 
absorbing configurations. 

4. Time scales and criticality in the movement of the random walker 

To find macroscopic physical quantities reflecting the time scales of the spectrum and 
the critical behavior at h c , we investigated the movement of the random walker starting 
from the specific initial positions U= (L — 1,1) (unfavored corner) or F= (L — 1,1) 
(favored corner). In the following we concentrate on the time dependence of the order 
parameter coordinate. Defining the average order parameter distance from the starting 
point as 



s(t) 



|< d(t) -d(0) >| 



(4.1) 



the start configuration corresponds to s = and the opposite corner of the triangle T L 
to s = 1. 



Start at unfavored corner 

First, we consider the random walker starting at the point U. We shall see that the 
movement from the unfavored to the favored corner is linked with the lower excitations 
of the hamiltonian. 
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s(t) 




Figure 7. Average distance s(t) from starting point U=(L — 1, 1). 

Data from simulation of the RWT model for a = h, h = 0.3 < h c . The solid curve 

represents the data for L = 100, the dashed for L = 600. 



Typical pictures of the walker's position s(t) are given in figures || 0> H an d ||- 111 
figure H we show the product L ■ s as a scaling function of t, in figure [7| the distance s as 
a scaling function of t/T short (both for /i below the critical point), and in figures |8] and ^| 
the distance s as a scaling function of t/L (at and above the critical point respectively). 
One observes the following properties of the scaling functions: 

• h < h c : 

The random walker stays near U for an average time exponentially increasing with 
L (Eq.( p.3|) ). For large L, the path near the unfavored corner (start) is described 
by the scaling function 

Ls(t) oc (1 - exp(-At)) . (4.2) 

Here we find A = 1.8(2) • rrii(h) for a = 0.25. Due to the distribution of flip times 
for finite L the walker leaves this regime at different times (figure |^) . 
For the approach to the favored corner, one observes a different scaling (figure 0). 
From the data, one finds that s is given by 

s(t, L) = goo - exp(-A — - ) , (4.3) 

short I, ) 

where 

/x . const. . d iS 

Soo = hm s(t) = 1 — + • • • . (4.4) 

In Eq.(|4.3|), the constant A = 1.15(3) is independent of a. The approach to the 
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s(t) 




Figure 8. Average distance s(t) from starting point U=(L — 1, 1). 

Data from simulation of the RWT model for a = j, h = h c . The solid curve represents 

the data for L = 100, the dashed for L = 400. 



favored corner scales with t/T shoTt . In the scaling of figure [7], the initial path 
(Eq.( |4.2| ) and figure |J) of the random walker cannot be seen because it is reduced 
to a point. 

h = h c : 

We consider the scaling regime t —>■ oo and L — > oo, taking t/L fixed. For 
t < T short /2 = L/2a, the distance from the starting point increases in leading 
order linearly with time. We find 

s(t, L ) = \ (it + $Vi + const. + • • • j , (4.5) 
where 5 = 0.17(2) independent of a. 

For large t > T s h or t/2, the random walker approaches s = 1 exponentially, as in the 
case h ^ h c . The scale, however, is the size of the system 

s(t, L) = Soo- (soo - a) exp(-K(a-^- - ^)) , (4.6) 

with k = 1.7(1) independent of a. Typical data is shown in figure |^. 
At the critical point, the initial path of the random walker Eq.( f4.5| ) is universal 
and linear in leading order. This reflects again the universal 1/L spectrum found 
for the critical point ( |2.6j ). It is remarkable that the exponential scale of Eq.( [4.6| ) 
is also universal (but not the scaling function). Up to the time normalization (as 
in section 2), they do not depend on a and scale with tL . 
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s(t) 




Figure 9. Average distance s(t) from starting point U=(L — 1, 1). 

Data from simulation of the RWT model for a = 4, h = | > h c . The solid curve 

represents the data for L = 100, the dashed for L = 400. The scaling variable is t/L. 



h > h c : 

Again we study the scaling t/L fixed taking t and L to infinity. Starting from 

the unfavored corner U, at first (for t < T s h or t) and in leading order the random 

walker follows a mean path determined by its velocity vector (see appendix A). The 

corrections are exponential in t. We find 

t C(l-exp(ft)) 
s(t, L) = ah- — . (4.7) 

With increasing h, the constant ( decreases and £ increases and thus the corrections 
become smaller, consistent with the explanation of appendix A. 
For t > T s hort, one again finds an exponential approach to the final value s = 1 (see 
figure |^) . One has 

X 

s(t,L) = -rjexp{-(mi + T)( t - T short)) , (4.8) 



where T short is given by (|3.4|). The scale is the gap m\ from the spectrum of the 
time evolution operator. Following appendix A, one finds the prefactor 
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Start at the favored corner 

We have also studied the path of the random walker starting from F. In that case the 
critical point h c does not play a special role. From the simulations, one finds 

s(t , L)oc i^PM, (4.10) 

where A = 1.25(5) • m 2 {h) for a = 0.25. This equation is similar to Eq.( |4.2j ), but note 
that in this case the movement is related to the time scale m 2 that describes a faster 
movement than the time scale m\. 

In summary, from the study of the movement of the random walker described in 
this section one can identify different scales. They are related to those excitations of H 
that describe the movement in question. One finds that the lower band of excitations 
(slow modes) is related to the path from F to U and the higher band (fast modes) to the 
movement around the favored corner F. This is confirmed by the results of appendix B 
for the model with absorbing configurations. For h c , we find a scaling function involving 
only the system size L as a scale, and depending on tL~ l . This scaling is analogous to 
the case of conformal invariance in critical equilibrium systems. 

We want to point out that, of course, the described scaling properties can also be 
found by investigating order parameters different from d. 



5. Conclusions 

In this paper, we studied the spectrum of the time evolution operator and other time- 



dependent properties of the two-dimensional non-local random walker (|1 . 1|) . This model 
is a simple example of a class of non-equilibrium models which exhibit spontaneous 
breaking of a symmetry. 

The reduced size of the configuration space of the single particle model of this paper 
allowed us to determine the spectrum numerically with enough precision to extract the 
gaps (time scales) and the large L limit of the low excitations. For many-particle models 
a numerical investigation of the spectrum is much more difficult and has not been done. 
The investigation of the spectrum of the RWT model led us to the following discoveries. 

Breaking the symmetry explicitly with a small h > 0, one is left with two phases, 
a meta-stable and a stable one, in accordance with the fact that the time evolution 
operator of the model has two ground states in the large L limit. Above the two ground 
states, one has two bands of excitations corresponding to two time scales. 

At the point h c = the meta-stable phase becomes unstable and, above that 
point, one is left with only one ground state and correspondingly one (stable) phase. 
This spinodal point has very interesting properties. Fixing h = h c , not only does 
one of the time scales vanish, but the corresponding excitations also acquire the form 
Ei = aki/L. Hence, it is universal in the sense that the constants ki do not depend on 
the remaining free parameter of the model. 
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At h — h CJ the dynamical critical exponent is z = 1, and the time scale vanishes with 
an exponent x = 2. Looking for "conformal towers", we did not find any regularities in 
the real part of the spectrum, but could show that the imaginary parts of the k t have 
equidistant seperations. We believe that the 1/L dependence on the system size is a 
more general phenomenon (z = 1 for the first excitation had already been found for 
the three states model in 0) and expect more structure to be found for the spectrum 
of many-particle-models at their spinodal points. The missing "conformal towers" for 
the RWT model may be a result of the reduced configuration space. Whether such a 
spectrum has further deep implications, as would be the case for equilibrium problems 
where one has conformal invariance, remains uncertain. 

Besides the numerical diagonalization, we performed simulations of the model to 
identify the physical relevance of the time scales. The lower gap is linked with the 
movement from unfavored to favored phase (slow modes), the higher gap with the 
movement within the favored phase (fast modes). We could identify different scaling 
regimes for the movement of the random walker starting at specific sites. The scaling 
functions reflect the time scales found for the hamiltonian. At the spinodal point, the 
appropriate scaling variable for the movement starting from the unfavored phase is t/L 
and the corresponding scaling function is universal. 

We also measured the first passage time for the random walker between stable 
and meta-stable phase. We found that, below the critical point, both these flip times 
increase exponentially with L. The flip times can be computed from the excitations 
of a modified model with absorbing configurations (appendix B). At the critical point, 
the first passage time from the unfavored to the favored corner of the triangle has a 
non-trivial, universal linear dependence on L, consistent with the dynamical critical 
exponent z = 1. The proportionality constant is different, by a factor of two, from its 
value above the spinodal point where the flip time also increases linearly with L. (There 
it can be calculated by a simple argument.) 

As with the three-state model, we have also investigated the free energy functional 
as a function of the order parameter d (|1.4j) . Surprisingly for a non-local model without 
detailed balance (but analogous to the three-state model), we observed a link between 
steady state properties (the free energy functional) and non-stationary dynamics (flip 
times). 

Those of the described results that relate to the region h < h c can also be found in a 
simplified, one-dimensional random walker model (RWRA model, see appendix C). For 
this model one can calculate the free energy functional, first excitations and flip times 
analytically. The further reduction of the configuration space leads to the different 
critical exponent x = 1 instead of x — 2. 
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Appendix A. Movement of the random walker for h > h c 

In the interior of the triangle 72, the random walker defined by ( |1 . 1| ) has an average 
velocity of 

" = ( ) (A.1) 

independent of the site (j,k). Here, the first coordinate is along the j-axis (number 
of — particles) and the second along the A;- axis (+ particles). The projection onto the 
direction of the order parameter d is given by 

v d = -2h(l-a) (A.2) 

Repeating the argument of Ref . , one sees that for small h both components of v 
are negative. The walker near either of the corners (L — 1, 1) or (1, L — 1) drifts towards 
the j- or /c-axis and subsequently jumps back to the respective corner. At h c = the 
velocity is antiparallel to the j-axis. For h > h c , the second component of v is positive, 
rendering the unfavored corner around (L — 1,1) unstable, because the random walker 
moves away from the j-axis and the non-local jumps along this axis become irrelevant in 
leading order. We have checked this with Monte Carlo simulations, (see also section ^). 

Neglecting these non-local jumps, from the mean path given by v one can easily 
calculate the time taken for a walker starting at U = (L — 1, 1) to hit the /c-axis and to 
reach the favored corner. One finds the flip time (|3.4j ). In the same way, (|A.2|) explains 
the leading term of (|4.7| ) and (|4.9|) . 



Appendix B. The RWT model with absorbing configurations 

In this appendix, we investigate further the relation of flip times and spectra. We focus 
on the case where the time that the walker remains in a phase increases exponentially 
with the system size. 

The flip times Ti ong (|3.2|) and T s h or t (|3.3|) are defined as first passage times. It 
is known that the corresponding time scales are related to dynamical properties, i.e. 
excitations, of a modified model with an absorbing state [[0|. An absorbing state is a 
configuration the system cannot escape from. Here we have to choose U (for Ti ong ) or F 
(for T short ) as the absorbing configuration. 



We shall first present the general argument (see [19|]). Consider a reaction diffusion 
system with an absorbing configuration {/? a b s } and described by a master equation ( |1.2j ) 
with a hamiltonian i7 a bs- The system may start in a configuration {fio} at time t = 0: 

\ P ^-°>W-{ for {(3} ^ {p } (R1) 

where the subscript of a vector denotes the component {/?}. Let E\ and \P\ > be 
eigenvalues and eigenvectors of the hamiltonian (Eq = being the ground state energy). 
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Note that the ground state |Po> is given by 

_ f 1 for {(3} = {/3 abs } 
l P o>{/3}-| for{/3}^{/W (R2) 

If i7 a bs is diagonalizable, there exist a A e C such that 

|P(t = 0)>=J> A |P A > . (B.3) 

A 

Obviously, the probability P in that the system has not reached its absorbing state 
is given by 

P in (i) = 1 - P({/WM) 

= 1 - (exp(-Pf abs t) \P(t = 0) >) {/W (B.4) 

= 1 -E°a exp(-P A t) |P A > {/W . 

A 

Denote as p(t)dt the probability that the system reaches {/5 a b s } during the time interval 
[t,t + dt]. Then 

p(t)dt = P jn (t)-P in (t + dt) . (B.5) 
Thus, is given by 

dP ft) 

pit) = — m = ~ £ exp (- EA *) ' Pa >w-> ' ( B - 6 ) 

A 

and its mean value, the flip time, is 

qa |P\>{/w 



T = —— iPabsl (B 7) 

A^O 



Pa 



From the last equation, it follows, under the conditions 

Pi oc exp(-m L) , lim P; > for i ^ 0, 1 , 3y > : lim P^ai 7^ , 

that in leading order 

T(L) oc exp(mL) . (B.8) 

The result is similar for any finite number of excited states Pj oc exp(— m^L). In that 
case, the scale of the flip time is given by the largest m^, with lirn^oo L%j 7^ for 
some y. 

On the other hand, if none of the excited levels approaches zero for L — > 00, one has 

T(L) < C (dim(L)) (B.9) 

independent of {fio} (dim(L) is the dimension of the configuration space for system size 
L). This implies that there is only one phase in the large L limit. 

We now turn to specific results for the RWT model. We have diagonalized the 
modified hamiltonian numerically. From the spectra we find the length scales /ii ong and 
/-ishort of the flip times: 
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• With the favored corner F made absorbing, the scale to ( |2.1| ) remains unchanged 
as does the gap of the lower band of excitations mi. This explains the relation 
/-^short = w ( P-7|) which is a special case of ( |B.8| ). Similar relations are found 
analytically for the simplified model discussed in appendix C. 

The lower gap vanishes at the critical point h c and the 1/L dependence of the levels 
is the same as for the non-absorbing RWT model (table p. The upper gap changes 
to a new value 1713 > m-i- 

These results confirm the interpretation that the second ground state (which exists 
for h < h c , L — > 00) and the lower band of excitations are linked with the meta- 
stable phase and with the movement from the unfavored to the favored corner. 
The corresponding scales of the spectrum remain unchanged, while the movement 
within the favored phase, and correspondingly the upper band of excitations, is 
altered by the absorbing configuration F. 

• In the case of an absorbing configuration U in the (originally) unfavored phase, there 
is more change to the spectrum. There is a first excited state E\ cx exp(— oj' L) 
with uj' 7^ to for any h, connected with two steady states for L — > 00. This reflects 
the fact that, with absorption in the corner U, both corners U and F constitute a 
phase for any h. From (|B.8| ) one has 

/•Hong = for any h, (B.10) 
analogous to ( |3.7|) . This relation is consistent with the numerical data to a precision 
of about 15%. 

Further, the scale uj' is related to vri2 of the RWT model without an absorbing 
configuration. For a = 0.25 one has 

uj' = 2.5(2) m 2 . (B.ll) 

As a further change to the spectrum, the lower band of excitations is shifted to a 
new band of excitations with a gap m 3 > m%, while the gap m-i remains unchanged. 
The new upper band of excitations is related to movement within the former 
unfavored corner. Note that m 3 (h) takes the same values with either the favored 
or the unfavored corner made absorbing. The corresponding eigenstates describe 
the (rapid) movement near an absorbing corner. 



Appendix C. Analytical results form the RWRA model 

In this appendix, we investigate the model defined by Eq.( |1.5| ). It is expected to give the 
physics of the RWT model below the spinodal point h c . We shall indeed see that for the 
RWRA model there are first passage times from one corner to the other which increase 
exponentially with L as for the RWT model. We calculate their exact expressions 
(Eqs. (|C.8|) and (|C.9|) ) and find the two scales with which these flip times diverge. 

We are also able to compute the first excitation (Eq.( |Cl5D ), which becomes a second 
ground state in the large L limit, and the first excitation for the RWRA model with 
either of the corners made absorbing (Eqs.( |CT2"0"D and ( |C21[ )). The first excitations are 
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proportional to exp(— ujL). In the case of no absorption or absorption at the favored 
corner, u is equal to the scale of the short flip time (|C.9|) ; and, in the case of the 
unfavored corner made absorbing, u is equal to the scale of the long flip time ( |C.9| ). 
All these results are analogous to those for the RWT model, but here they are obtained 
analytically. 

We have also computed the free energy functional Eq.( p.26 ), whose graph is very 



similar to that for the RWT model's. As in the RWT model, the height differences 
/max — /min give the scales with which the flip times (|C.8HC.9|) diverge. 

Appendix C.l. Flip times 

As for the RWT model, for h = the random walker stays most of the time near the 
points F=(l, L — 1) and U=(L — 1, 1). For h > 0, the point F is favored and the stable 
phase is concentrated around this site. 

To calculate the flip time from the meta-stable to the stable phase, we define Dj^ 
as the first passage time from (j,k) to (1, L — 1). 

Concentrating on the first step of the random walker (in an infinitesimal time interval 
dt), we find 

D lifc =(l-dt)p a , fc + dt) 

+ ^-dt(D 1>k _ ± + dt) 

+ I±^dt(£> liL _i + dt). 
for j = 1, k = 2, . . . , L - 2 (cf. §). This yields 

D l>k = 1 + ^£>i, fe -i + (C.l) 
Similarly, the other conditions for the Dj^ read 

Dj,i = 1 + ^^-1,1 + j = 2, . . ■ , L - 2 (C.2) 

Di,x = 1 + ^D l>L _, + ^Dl-i,i (C.3) 

1- h 1-h N 

— - — D\,L-1 = 1 + -^-A,l-2 (C.4) 

From the definition of Dj^, the boundary conditions are given by 

£l,L-l = 

Dl-i,i = T short . (C.6) 



The recurrence relation (|C.2|) can be solved as 



2 / [1 + hV 2 fl + h^ j 



d ^ = d "—h[{— - — +1 • (C7) 
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and the condition ( |C.5| ) gives us the flip time 

r ~ = (t^)' (its)" - ■ < c -«) 

Replacing h by — h, we obtain the flip time from favored to unfavored phase 

One sees that the flip times diverge exponentially with L. The scales are /z s hort = 
log(2/(l — h)) and /ii ong = log(2/(l + h)). This implies that there are two ground 
states in the large L limit. Linear combinations of these give the phases; a stable phase 
concentrated around F= (1, L — 1), and a meta-stable around U=(L — 1, 1). 

Appendix C.2. L dependence of the first excitation 

With the notation x s = P(s, 1), y s = P(l, s) and x% = y\ = P(l, 1), the master equation 
of the RWRA model reads 

1 i u 

s = 2,...,L-2 (CIO) 
s = 2,...,L-2 (C.ll) 



d 


-X s + 


1 + h 


~&t xs ~ 


2 Xs+1 


d 




1 - h 


dt ys = 


Vs + 


2 Vs+i 


d 




1 + h 


dt xl = 


—Xi + 


Xo + 

2 


d 




1 - h 




= -XL 


1 + 2 








d 




1 + h 


dt^- 1 


= -Vl- 


- 1+ 2 



V^ 2 (C12) 

X> ( c - 13 ) 

i=l 

eV ( c - 14 ) 

8=1 

Setting x s = exp(— Et)X s and y s = exp(— Et)Y s , and solving the resulting recurrence 
relations, we find the condition on an eigenvalue E, 

o = e(e(i- e) 2L - 3 + (1 ~ L ft _ 2 1 )I '' 



2 L 



,2L-sA + h ^A-h 



-1 



^(1 - EY L - b {—^- - E){— E)j . (C.15) 

There is one solution with E = 0. The corresponding eigenvector is the steady state 
for finite L. For L ^> 1, the first excited level behaves as 

E x = \{l - h 2 ) (xexp(^L). (C.16) 

Thus, for L — > oo, one has two stationary states resulting in two phases. 



23 



The higher spectrum has a two band structure similar to that for the RWT model. 
The corresponding gaps are mi = (1 — h)/2 and m 2 = (1 + h)/2. Thus, one finds 
uj = rrii for h — > 1. (For the RWT model one has proportionality of uo and mi for any 



h, see ( |2.4j) .) We point out that for the RWRA model the gap vanishes with a critical 
exponent x = 1, while for the RWT model one finds x = 2 (see Eq. fl2.4| )). This is due 
to the reduction of the configuration space. 

One can compute the large L behavior of the flip times from the L-dependence of 
the first excited level of a modified model (appendix B). For the flip time T s hort, one 
implements an absorbing site at the stable corner (1, L — 1), replacing equations ( |C.11| ) 
and qCTHD by 

yVu s = 2,..., £-3 (C.17) 

(C.18) 

eV ( C - 19 ) 

With this, one finds a first excited level 

E (l,L-l)abs. = 1 (1 _^ 2) ^1+^^ L \ (C20) 

With (L — 1, 1) absorbing, instead of (1, L — 1), one analogously obtains 

£(L-l,l)abs. = 1 (1 _ h2) L 2 . (C.21) 



d 


1 


dt vs = 


Va + 


d 




dt yL - 2 


= -J/1,-2 


d 


1 + h 


d^- 1 


2 



This way, with equations (|C.8|) , (|C^9|) and ( B.8|) , we find 

2 N L ~ 2 

-^short 0^ 



2 ^ L " 2 



.1 - 

which is indeed consistent with the flip times ( |C.8|) and (|C.9|) computed directly above. 
Appendix C.3. The stationary probability distribution 

To calculate the probability distribution of the steady state, one has to solve equations 
(PTOD -(|OT4|) for f t x s = ±y 8 = 0. One finds 

X s = \n{1 - h 2 ) f-^-V (C.22) 



4 v ' \ l + h 



Y s = \N(1 - h 2 ) (— S = 2,...,L-1 (C.23) 
X 1 = Y 1 = N, (C.24) 
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where the normalization N is 

(2(1 - h)) L ~ 2 + (2(1 + 1W' 2 - (1 - h 2 ) 
The free energy, defined by (|1.3|), takes the form 



/(d) = < ~ (C26) 

log^+dlog^) forrf<0. 

The minima of the free energy functional are at d — — 1 and c? = 1, corresponding to the 
stable and metastable phase respectively. In between, the free energy functional has a 
constant slope. 

In conclusion, we have shown that the RWRA model presented in this appendix 
captures most of the physics of the RWT model below the critical point. 
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